Reply on "Fluctuation-dissipation considerations for phenomenological damping 
models for ferromagnetic thin films" [N. Smith, J. Appl. Phys. 92, 3877 (2002)] 
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We show that the critique of our recent papers presented 
in the abovementioned paper (NS) appeals to an incorrect 
mathematical analogy between electrical circuits and linear 
magnetization dynamics, improperly uses classical concepts of 
normal modes and basic equations, gives inconsistent results 
and therefore comes to incorrect conclusions. 

The study of linear stochastic magnetization dynamics 
is of great importance in applications to nano-magnetic 
devices based on ultra-thin films. Gyromagnetic mag- 
netization motion around an effective field is randomly 
forced by fluctuations on spins by means of interaction 
with a thermal bath (phonons, magnons, conduction elec- 
trons, impurities, etc.). The role of the thermal bath can 
be approximately reduced to terms describing relaxation 
and random fields in the magnetization dynamic equa- 
tions. 

For simplicity, the conventional theoretical approach 
is purely phenomenological. ft is based on the Landau- 
Lifshitz equation [|l| with random fields, or its modifi- 
cation in a Gilbert form This approach has the de- 
fect that the phenomenological damping term in Landau- 
Lifshitz (and Gilbert) equation was introduced just for 
the case of high magnetic symmetry (axial symmetry). 
An applicability of an isotropic damping together with 
corresponding noise to strongly anisotropic system such 
as thin film has not been proven |^]. 

In our recent papers Q, |^ we have developed a the- 
oretical approach based upon the representation of the 
magnetization dynamics as the motion of damped non- 
linear oscillator driven by a random force (thermal fluc- 
tuations). The oscillator model is a convenient tool H] 
to establish a "bridge" between the microscopic physics, 
where the oscillator variables a* and a naturally describe 
spin excitations (as creation and annihilation operators), 
and the macroscopic stochastic magnetization dynamics 
for normal modes. We have calculated the magnetization 
noise spectrum in a thin film for a physical loss mech- 
anisms [0 and shown that our result differs from that 
of obtained by the phenomenological approach and is in 
agreement with experiment 

Recently Neil Smith (the paper cited in the title, later, 
NS) criticized our approach claiming that "a proper un- 
derstanding of the FDT (fiuctuation-dissipation theo- 
rem) ... may have been complicated by several recent 
papers Q], l^". In this Reply we argue that this critique 
1) appeals to an incorrect mathematical analogy between 



electrical circuits and linear magnetization dynamics, 2) 
improperly uses classical concepts of normal modes and 
basic equations claiming simultaneous diagonalization of 
three hermitian matrices, and 3) gives inconsistent re- 
sults. 

1. Incorrect mathematical analogy. Any analogy in 
physics is based on a similar structure of mathematical 
equations. The basic equations for linear electric circuits 
have the form (here and later we use notations of NS): 
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This is a set of the second order differential equations, 
where the scalar variables Xj describe electric charges, 
Fi{t) are external voltages and symmetric matrices de- 
scribe inductances (Mij), resistors (Dij) and capacitors 

On the other hand, linear magnetization dynamics do 
not contain inertial terms {Mij = 0) and do contain an 
antisymmetric gyromagnetic matrix Gy : 
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This is a set of first order differential equations and there 
is no mathematical analog for Gij in a linear electric 
scheme. In order to describe N coupled magnetic oscil- 
lators (spin waves) each dynamic variable Xj should be 
chosen as a column (not a scalar), containing two trans- 
verse magnetization components in the j-th micromag- 
netic cell: 
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In this formulation there is no direct mathematical 
analogy between linear electrical circuits (0) and linear 
magnetization dynamics (|^). Thus, without loss of gen- 
erality, the compilation of "electrical analogies" (Sec. Ill 
in NS) has no mathematical consequences for magnetic 
systems 0. 

2. Incorrect claim of simultaneous diagonalization of 
three hermitian matrices (including damping D ). We 
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think that this claim is just a simple misunderstand- 
ing in NS of the normal mode (spin wave) concept in 
magnetic dynamics. Spin waves are a principal concept 
of a conservative magnetic system where the energy of 
the linear magnetic oscillations can be diagonalized and 
represented as a set of independent harmonic oscillators 
(see, e.g., Uniform rotation of the magnetization 

corresponds to a spin wave with zero momentum. Spin 
waves interact with each other and with other degrees of 
freedom (elastic waves, conduction electrons, etc.). Spin 
wave damping appears as a result of reducing (averag- 
ing) these interactions (interaction with a thermal bath). 
Therefore the reduced equations can not have a general 
structure. The damping matrix is specific to the mag- 
netic system and can not have an arbitrary form, as as- 
sumed in NS. 

3. Inconsistency of obtained "general" results. Let 
us check the general results obtained in NS for the 
fluctuation-dissipation theorem as used in Gilbert and 
Bloch-Bloembergen dynamics. In NS the Gilbert equa- 
tion is written in the form: 
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On the other hand, in NS the Bloch-Bloembergen 
equation has been written as: 

- G* + (*H - G* /T2) ■ = h{t) (10) 
dt 



with corresponding FDT of the form: 
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We see that Eqs.(|9|) and (|10| ) have similar forms. They 
coincide, for example, in the case 
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and 1/72 = a-fHo. In this case their FDT relations (g) 
and (|ll| ) are distinctly different from each other. This 
fact simply indicates that NS results are inconsistent. 

This work was partly supported by matching funds 
from the Center for Magnetic Recording Research at the 
University of California - San Diego and CMRR incorpo- 
rated sponsor accounts. 



and H describes effective magnetic fields. The corre- 
sponding thermal fluctuation fields are given by: 
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referred as FDT in NS. It is well known that the Gilbert 
equation can be rewritten in a mathematically equivalent 
form of the Landau-Lifshitz equation. This procedure is 
most simple in the case of small damping when we can 
neglect the higher order damping terms (such as a^). We 
can find the magnetization derivative on time from Eq. (|^) 
without damping 
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and represent the damping term as 
dm 
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G j =7 G . Neglecting small a^f G • h (i), 
the Eq.(^ may be written as: 
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